Structure of unstable light nuclei 



D.J. Millener^ 

Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA 



o . 

O . Abstract 
(N 

^ . The structure of light nuclei out to the drip lines and beyond up to Z ~ 8 is 

interpreted in terms of the shell model. Special emphasis is given to the underlying 
supermultiplet symmetry of the p-shell nuclei which form cores for neutrons and 

■ protons added in sd-shell orbits. Detailed results are given on the wave functions, 
widths, and Coulomb energy shifts for a wide range of non- normal parity states in 

> ■ the p-sheU. 
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Ch ' 1 Introduction 

> 

X 

' The light nuclei have long provided a testing ground for nuclear models. Most 

^ of the basic ideas were already in place by the end of the 1950 's as can be 

seen from a perusal of the proceedings of the Kingston conference in 1960 [1] 
and a recent history of the development of our understanding of nuclear struc- 
ture by Wilkinson [2] (this reference lists the important papers and outlines 
subsequent developments). 

As far as the shell model was concerned [3], the fractional parentage coefficients 
for the p shell in both jj and LS coupling had been available since the early 
50 's but the slow development of electronic computing meant that it was the 
mid 50's before the first diagonalizations (up to 22 x 22) for {sdY by Elliott 
and Flowers [4] and for the complete p shell by Kurath [5] were published. 
Excitations across major shells with proper elimination of spurious centre-of- 
mass states were considered by Elliott and Flowers for A = IQ [6]. This paper 
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also provided a microscopic description of giant dipole strength in ^^O. The 
presence of non-normal-parity states at low excitation energy was known in 
^^F at 110 keV (suggested as the weak coupling of a Opi/2 proton hole to ^°Ne 
[7]), in ^^N at 2.4 MeV (the weak coupling of an lsi/2 proton to ^^C [8]), and 
a I"*" ground state for ^^Be was strongly suspected [9,10]. In the appropriate 
limits, the close connection between the shell model and the cluster model and 
between the shell model and the collective (Nilsson) model was known, the 
latter being made clear from Elliott's formulation of the SU(3) shell model 
[11]; also the need for and basic origin of effective interactions and charges. 

Of course, the information on spins, parities, and other properties of nuclear 
levels was still rather sparse [12]. The next few decades brought a rapid filling 
out of this basic knowledge, as may be traced through the succeeding tabula- 
tions [13]. In recent years, the advent of radioactive beam facilities has greatly 
improved our knowledge of the structure of nuclei far from stability, out to 
the drip lines and beyond in the case of light nuclei. 

The purpose of this paper is to interpret the structure of light nuclei up to 
Z ~ 8 in terms of the shell model, using where necessary its relation to the 
cluster, Nilsson, and weak-coupling models to provide a simple understanding 
of the structure. In principle the shell model is complete and in practice it 
can be used phenomenologically to correlate rather precisely the systematics 
of nuclear properties over considerable ranges of mass number. For some mass 
numbers, the tabulations are up to 12 years out of date but the general flavour 
of advances in the field can be obtained from recent major conferences such as 
ENAM98 [14], recent reviews on specific topics [15-18], and the most recent 
references to specific nuclei (given later). 



2 Structure Calculations 

Large scale Ohw shell-model calculations have now reached close to the middle 
of the pf shell [19]. These results, obtained using a G-matrix interaction with 
monopole modifications, have also been used to benchmark shell-model Monte 
Carlo calculations [20,21] which can be used for even larger model spaces. Both 
types of calculations have been employed for neutron rich nuclei involving both 
the sd and pf shells [20,22] where some ground states involve the excitation 
of pairs of neutrons across the A'" = 20 shell closure. Similar violations of the 
normal shell ordering occur at = 8 [23] . As is discussed later, the structure 
of such nuclei can be described by shell-model calculations which use two-body 
matrix elements fitted to O/lo;, Ihuj, and 2fuj configurations in the A — 10 — 22 
nuclei [24]. 

Binding energy effects are not taken account in the shell-model calculations 
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themselves. Rather, Woods-Saxon radial wave functions evaluated at the phys- 
ical separation energy are used to calculate transition matrix elements. Further 
improvement is possible by using the shell-model structure information to set 
up radial equations for the appropriate one-nucleon overlap functions [25] . 

Microscopic cluster models generally include the correct degrees of freedom to 
describe extended, loosely bound systems without violating the Pauli princi- 
ple or introducing spurious center of mass excitations. Recent developments 
include solutions by the stochastic variational method [26], antisymmetrized 
molecular dynamics [27], molecular orbital methods [28], and generator cor- 
dinate methods [29]. Such calculations use saturating central forces, such as 
the Volkov or Minnesota interactions, augmented by a spin-orbit interaction. 
Parameters, such as the space-exchange mixture, are often varied on a case- 
by-case basis to reproduce energies with respect to thresholds for states of 
interest. 

Finally, there have been great advances in the theoretical treatment of few- 
nucleon systems [30] and some of the techniques can be used for p-shell nuclei. 
In particular, variational Monte Carlo and Green's function Monte Carlo re- 
sults using realistic free NN interactions and phenomenological but theoreti- 
cally motivated NNN interactions have been published up to A = 8 [31]. Very 
impressive agreement with experiment has been obtained for ^Li(e, e')^Li form 
factors [32] and for ''Li(e, a'p^Re momentum distributions and spectroscopic 
factors [33]. The last result emphasises the essential quasi-particle nature of 
the shell model and the role played by (short-range) correlations [34,35]. 



3 The p-shell nuclei 

A comprehensive description and understanding of the structure of p-shell 
(Ohw) states is important both in its own right and because neutrons in sd- 
shell orbits are added to p-shell cores as one moves towards the drip line or 
higher in excitation energy for many light elements. In fact, there is often near 
degeneracy, or coexistence, of nominally Ohw, Ihw, 2hu, etc. states. 

It should be noted that Barrett and collaborators have performed ab initio no- 
core shell-model calculations in at least 4/ia; spaces up to A = 12 with effective 
interactions derived microscopically from realistic NN interactions [36]. The 
spectra for known p-shell levels are good. To date, however, the non-normal- 
parity states, and consequently all multi-/ia; states, are predicted too high in 
energy. 

At the beginning of the p-shell, the wave functions for the observed states 
have long been known to be close to the supermultiplet (LS) limit in which 
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the orbital wave functions are classified by the SU(3)dO(3) quantum numbers 
(A ji)KLL (A = /i — /2 and /i = /2 — /s from the spatial symmetry [/] = 
[/1/2/3]) and the SU(4)dSU(2) xSU(2) quantum numbers [J]I3TS (see, e.g., 
Barker's wave functions for A = 6 — 9 [37]). In fact, the supermultiplet scheme 
provides an excellent basis for an understanding of the structure and energetics 
of all p-shell nuclei. This can be demonstrated from an analysis of a p-shell 
Hamiltonian obtained by fitting 34 levels for A = 10 — 12. In the fit only 
a limited number of well-determined linear combinations of the parameters 
were allowed to vary. The single-particle energies were always well determined, 
as were the central matrix elements with the exception of the singlet-odd 
interaction. The strength of the tensor interaction was fixed to obtain, in 
competition with the spin-orbit interaction, the sign and magnitude of the 
and ^Di mixing in the ®Li and ^"^N ground states necessary to explain 
the small Gamow-Teller matrix element in ^^C (3 decay and a small negative 
quadrupole moment for ^Li. The spin-orbit sphtting of 4.78 MeV (3.49 MeV 
for a similar A = 6 — 9 fit) is much larger than that of the Cohen-Kurath 
interactions [38] for the light p-shell nuclei. 

The calculations were performed with an SU(3) shell-model code based on 
the formalism of French [39] adapted for SU(3) [40]. Prom the coefficients 
of the six (0 0) and (2 2) SU(3) tensors which define the central interaction, 
it is possible find the coefficients of the 5 SU(4) invariants which character- 
ize the symmetry-preserving part of the interaction [41]. This decomposition, 
including the contribution from the one-body centroid, is given below. The 
coefficient of the sixth tensor is essentially zero so that the central interaction 
is almost perfectly symmetry preserving. 

H = 1.563n - 1.709/ij - 3.908Pij- + 0.590L2 - 1.0785^ ^ 0.594X2 ^^-^-^ 

For the spatial symmetry [/] = [/1/2/3], the expectation value of the space- 
exchange interaction Pij is | I]j /«(/« — 2i + 1) = 11^ — Ua, where Ug — ria is the 
difference between the number of symmetric and antisymmetric pairs. 

The binding energies given by Eq. (1) for the lowest p-shell states in nuclei 
of interest are shown in Fig. 1. As can be seen, the dominant contributor to 

the spacing between the nuclei with different isospins is the space-exchange 
interaction. These spacings are more or less conserved when the non-central 
interactions are switched on with admixings of < 30% (often much less) of 
lower symmetry into the lowest states throughout the p shell. For example, 
the ground state of -^^C is typically ~ 80% [44] symmetry (LS = 00), as we 
might expect for a nucleus that we can think of as 3 a particles, with a 
[431] admixture (LS = 11) of ~ 16% induced by the spin-orbit interaction. 
This admixture provides the Gamow-Teller (GT) transition strengths in the 
P decay of ^^B or ^^N. The closed P3/2 shell forms < 40% of the ground-state 
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Fig. 1. The solid lines labelled by spatial symmetry, nucleus, and the spin-parity of 
the lowest p-shell state give the Coulomb-corrected binding energies with respect 
to "^He for the central part of a p-shell interaction fitted to 34 level energies from 
the ^ = 10 — 12 nuclei. The dotted lines represent the corresponding experimental 
binding energies for the lowest state of each mass number. These energies are well 
reproduced when the non central interactions are included, with the dominant effect 
(typically an energy gain of ~ 4 MeV) coming from the one-body spin-orbit interac- 
tion. The circled numbers give differences in the eigenvalues of the space-exchange 
operator for the spatial symmetries [/] shown. 
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wave function (the overlap of p|/2 [44] is only y5/81 [8]). 



With regard to the smaller angular momentum and isospin dependent terms 
in Eq. (1), the term gives the well-known L{L + 1) dependence exhibited by 
the p-shell interaction. The influence of the and terms can be seen in the 
strong low-cncrgy Gamow- Teller transitions observed in the high Q/j decays 
of ^He, ^Li, and ^^Li. In the case of ®He decay {Qp = 10.653 MeV) [42], the 
lowest three 1"*" states of ^Li have small GT matrix elements because the wave 
functions have dominantly [31] spatial symmetry (mixed LS = 10, 11, 21) but 
the I4 state near 9 MeV has a large GT matrix element corresponding to the 
transition I [22] L=0S=0T=2) |[22] L=QS=1T=1). The intensities of these 
configurations in the initial and final states are 74.0% and 71.5%, respectively, 
and the influence of the and terms in lowering the I4" state below the 
analog of the initial state at 10.822 MeV in '^Li is clear. Similarly, in the decay 
of ^■'^Li {Qp — 20.7 MeV) [43,44], the strong transition is between states of 
[322] symmetry with L = 1 to a state just above 18 MeV in ^^Be. The main 
candidate is a | state predicted at 17.8 MeV with 87% by intensity of [322] 
symmetry with S* = |. An essentially degenerate | level is also predicted to 
carry significant GT strength. The predicted p-shell GT values are strongly 
quenched by large (so?)^ admixtures in the ^^Li ground state [45,46]. 

Another important aspect of the underlying supermultiplet symmetry for p- 
shell states is how it influences parentage, single-nucleon or multi-nucleon, 
and therefore the magnitudes of cross sections for transfer reactions and the 
widths of unbound levels. Before turning to p-shell properties of interest for 
some of the He, Li, and Be istopes, a brief discussion of these matters is given 
in the following subsection. 

3.1 Parentages and widths 

For unbound states with modest widths an excellent way to estimate the 
nucleon decay width at energy above threshold is given by 

TM^C^Sr^iEN) , (2) 

where Ej^i — i^V^^ is the complex resonance energy for a potential well [47], e.g. 
a Woods-Saxon well with vq ~ 1.25 fm and a ~ 0.6 fm, with a depth chosen to 
reproduce E^. The spectroscopic factor S is the square of the matrix element 
of a creation operator between initial and final states and C is the isospin 
Clebsch-Gordan coefficient. In a schematic form, between basis states for n 
particles in a shell, 

{n\\\a'^\\\n — 1) = ^/n{n — l,l\}n) (3) 
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where the coefficients of fractional parentage (cfp's) are products of weight fac- 
tors and Clebsch-Gordan coefficients for SU(3)dO(3) and SU(4)dSU(2) xSU(2) 

[48] . The weight factors ■^j'^ are particularly important because they give the 

branching to states of final symmetry [/'] which can be reached from an ini- 
tial symmetry [/]. If more than one [/'] is allowed, there will be parentage to 
widely separated states of the A — 1 core nucleus. For example, for A = 13 
the weights for [441]— i^[44] and [431] are and y^, respectively. On the 

other hand, nucleon removal from ^^C leads only to the f^^ , |^ , and 1^ states 
of ^■'^C or ^^B because [44]— >[43] is the only possibility. For N > Z, the sum 
rule S{Ty) = Z applies for proton removal, with the T> states generally 

belonging to a different [/'] from the T< states, while both T< and T> states 
contribute for neutron removal. 

It is very easy to generate spectroscopic amplitudes for multinucleon transfer 
from the one-particle cfp's of Eq. (3). For 3 and 4 nucleon removal, it is 
assumed that the k nucleons have maximum spatial symmetry [k], which can 
be projected onto a (Os)^' internal wave function for the cluster. Amplitudes 
for 1 — 4 nucleon transfer have been tabulated [49] for the Cohen and Kurath 
wave functions. The generalization to other major shells and mixed major 
shells is straightforward. 



3.2 The He isotopes 



There is now considerable evidence for a number of relatively narrow excited 
excited states in ^~^°He [18]. The positions of the lowest narrow states in ^He 
and ^°He are consistent with the p-shell expectations of Ref. [50] and this work. 
The widths of ^He(|^), ^He(|^), and ^He(2|) calculated according to Eq. (2) 
are in good agreement with the experimental values. However, the nominally 
Pi/2 state {S{th) = 0.65 - 0.83) at 1.13 - 1.27 MeV above threshold in ^He 
[18] is predicted to have a width in excess of 1 MeV, which more than a factor 
of 2 larger than is reported. Very recently, strong evidence has been obtained 
that the ground state of ^He is an unbound s-state with a scattering length of 
< -10 fm [51]. The WBP and WET interactions [24] predict that the state 
should be lowest [51] followed by a | state, and a ^ ) § pair 2.2 and 2.4 
MeV above the |^ state; Poppelier et al. make very similar predictions [52]. 
The 1^ state could correspond to the second state seen in the ^Be(^^C,^^0)^He 
reaction [18] at = -S^ ~ 2.4 MeV (Tf = 624 keV). In any case, transfer 
reactions of this type should populate a number of narrow states of the form 
^}ie{0~^ ,2~^)0{sd). The presence of the low-lying sd states in ^He clearly has 
implications for the structure of ^°He. 

Recently, an excited state of ^He which decays mainly into the ^He-|-3n channel 
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has been observed via the p(^He, c^)^He reaction [53] at = 2.9(3) MeV with 
r = 2.2(3) MeV. The p-shell | state, which is predicted shghtly higher in 
energy, has the correct decay properties but must be produced via a two-step 
process [53]. Also, a broad state at ~ 3.2 MeV is seen in the '^Be("^^N,^'^F)'''He 
reaction and assigned — ^ [54]. However, a | state at this energy would 
be exceedingly broad and, in addition, the calculated spectroscopic factor for 
two-proton pickup from ^Be is vanishingly small. The | state, on the other 
hand, has a substantial L=2, S—0 two-proton spectroscopic factor. 

3.3 The Li isotopes 

There have been many investigations of the unbound nucleus ^'^Li. These are 
summarized in recent works which confirm the existence of s-wave strength 
near threshold [55] and a p-wave resonance at 500 keV above threshold with a 
width of 400(60) keV [56]. Shell-model calculations predict a low-lying l"*", 2"*" 
doublet separated by only 23 keV for the A = 10 — 12 interaction (and 169 keV 
for the A = 6—9 interaction) with S = 0.98 and 0.71, respectively. At En = 500 
keV, the predicted width for the 1+ state is 375 keV; a shghtly narrower, 
nearly degenerate 2'^ state could also be populated in the experiments. The 
^•^80(^^0, ^^N)^°Li reaction shows evidence for states at E^ = 0.24, 1.40, and 
4.19 MeV [18]. 

Based on the evidence from ^'^Li and general systematics, states of the form 
^Li{Ofiw)i^{sdY and ^^Li{Ohuj)0{sd) are expected close in energy to the simple 
p-shell I configuration for ^^Li. Indeed, the interactions of Warburton and 
Brown [24] predict essential degeneracy of the lowest Ohuj and 2huj configura- 
tions in ^^Li (and ^^Be). The evidence from ^^Li (3 decay for strong mixing has 
already been mentioned and an analysis of a fragmentation experiment [57] 
comes to the same conclusion. There is also evidence for a number of relatively 
narrow excited states [18] which are necessarily not of p-shell origin and must 
be explained by any realistic structure calculation. 

3.4 The Be isotopes 

Many years ago, Kurath and Pieman [58] showed that the lowest p-shell states 
generated using a central plus one-body spin-orbit interaction had an essen- 
tially perfect overlap with states projected from a Slater determinant of the 
lowest Nilsson orbits for a Nilsson Hamiltonian with the same spin-orbit in- 
teraction as the shell model and a deformation that varied in a very regular 
way throughout the p-shell. In cases where two bands with different K lie 
close together, linear combinations of the two J-projected states reproduce 
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the two corresponding shell-model states. In Elliott's SU(3) model [11], states 
with good orbital angular momentum are obtained by projection from a Slater 
determinant, or linear combination thereof, of asymptotic Nilsson orbits. In 
the case of non-zero intrinsic spin, Elliott and Wilsdon [59] projected states 
with good J and Kj — Kl + Ks from a product of the SU(3) intrinsic state 
and the intrinsic spin wave function. The J-projected states correspond to a 
well-defined linear combination of L-projected LS-coupling states. In the shell 
model, the spin-orbit interaction drives the mixture of states with different L. 
Of course, the spin-orbit interaction also mixes different SU(3) representations 
through its (1 1) tensor character. In the supermultiplet, or SU(3), basis it is 
very easy to see the band structure and to check the degree to which Kj is a 
good quantum number. 

The first K — ^ orbit for prolate deformation is filled at ^Be and the interplay 
of the remaining K — ^ and K — ^ orbits can be seen beyond A = 8. ^Be 

has a. K — ^ ground-state band with L — 1 and L = 2 in the ratio to 

— in the SU(3) limit (the same is true for -"^^B). For A — 10, two particles 
in the K = ^ orbit give rise to the K = 3 ground-state band of for T = 
(S* = 1) and to the K = ground-state band of ^'^Be for T = 1 (S* = 0). In the 
latter case, promoting one neutron to the next K = ^ orbit gives the K = 2 
bandhead at 5.958 MeV, as shown in Fig. 2. The reason for the very strong 
population of this state via a Gamow- Teller transition in the (t,^He) reaction 
on ^°B [60] can be seen from the L decomposition of the pure Kl — 2, Ks — 1 
^°B ground state, 

|'»B 3+ if. = 3) = = 2) - = 3) + = 4) , 



as can the reason for the excitation of the K^ = 2, 3+ state at 9.40 MeV. 
The 3.37-MeV 2+ state is seen weakly, some K mixing of the 2+ states being 
necessary to ensure the near equality of the 2f 0^ E2 (isoscalar) transitions 
in ^°Be and ^''C [13]. The 2+ state near 9.6 MeV in Fig. 2 is seen in proton 
pickup from ^^B [13] and corresponds to a p-shell state with [33] symmetry. 
States with [411] symmetry are also predicted near 10 MeV but they would 
have very large neutron decay widths. The near degeneracy of the Ihuj and 
2huj configurations has long been known [13] and is nicely reproduced in a 
number of recent cluster model calculations [26-28] . 

For the last three neutrons in ^^Be, one expects bands formed from {K = |)^ 
(K = |) and {K = ^){K = intrinsic states. In the shell-model the corre- 
sponding bands in the SU(3) limit have (A/i) = (2 1) ([421] spatial symmetry) 
with Kl = 1 and Ks = =t|; the details are shown in Table 1, where it can be 
seen that the wave functions are indeed dominantly of [421] symmetry with 
K — ^. With respect to identifying experimental candidates for these states. 
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Fig. 2. States of ^°Be below 10 MeV excitation energy. 
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Fig. 3. States of ^^Be below 4 MeV excitation energy. 
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Table 1 

States in the p-shell K = ^ band of ^^Be with the bandhead at 0.320 MeV for the 
fitted A = 10 — 12 interaction. 
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the tabulation [61] is misleading since the spin and parity assignments are 
based on {t,p) angular distributions contradicted by later work [62]. Pickup 
reactions from heavier p-shell nuclei provide a clean way to identify domi- 
nantly p-shell states. In the ^^C(^Li,^B)^^Be reaction [63], the | state and 
states at 2.69 McV and 4.0 MeV arc seen. Also, four T = | states in ^^B have 
been identified via the ^^C(p, a)^^B reaction [64], the first three of which were 
identified with first | , | , and | states, the latter two at 2.37 MeV and 
3.58 MeV relative to the first. The shell-model spectroscopic amplitudes are 
consistent with the assignments from these two reactions. The strong branch 
to 2.69-MeV level in the (3 decay of ^^Li [44] is model-independent evidence 
of negative parity for this level. 

The level scheme for states below 4 MeV shown in Fig. 3 is the same as that 
from Liu and Fortune's analysis of the ^Be(t,p)^^Be reaction [62] with the 
exception of the 3.956-MeV level. However, the L — 2 angular distribution 
observed for this level is consistent with db "2 clSSl gnment and the discussion 
in the previous paragraph indicates that one of the two levels in this region is 
the first p-shell | level. Charge-exchange reactions also see strength in this 
region [60]. The ratios of strengths to the 0.32 and 2.69 levels, and to the 3.9 
MeV region in ^^Li(/3~) decay [44] is consistent with the shell-model prediction 
although there is some uncertainty over the the division of strength between 
the 3.89 and 3.96 MeV levels. The authors of Ref. [44] favour an interchange of 
the I and |^ assignments for the 3.41 and 3.89 MeV levels in Fig. 3. However, 
the width of the 3.41 MeV level is over 100 keV [62] and this is inconsistent 
with the small ^3/2 spectroscopic predicted for the lowest |^ level and the 
calculated d-wave single-particle width of 144 keV. The predicted width for 
the p-shcll I level at 3.96 MeV is 17 keV for S(2+) = 0.66 compared with 
the experimental value of 15(5) keV; a slightly smaller value S(2+) = 0.49 is 
required to give the observed width of 201(10) keV for the apparent analog in 
i^B. 

All three states of the K'^ = | band have large parentages for neutron removal 
from ^^N. The proton spectrum which results from the decay of the resulting 
^^N states has been studied [65] and can be qualitatively reproduced using the 
calculated decay widths for the p-shell states of ^^N. 
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An important point is that the parentage for neutron removal from the [421] 
symmetry p-shell states is complex because the parentage to the allowed sym- 
metries [42], [411], and [321] is divided in the ratio 9:10:16. The model inde- 
pendent sum rules for S(r<) and S(T>) are 4.5 and 2.5. Of the C^S = 5 for 
neutrons, 0.5 must go to T = 2 states at ~ 22 MeV excittion energy in ■'^'^Be, 
1.8 to [42] states, and 2.0 to [411] states. For the lowest | state, about 2.2 of 
the parentage goes to the ground state and the first two 2+ states and about 
1.6 to states in the 9 — 13 MeV energy range, so that there is a small redistribu- 
tion of strength from the pure symmetry limit. This parentage is necessary for 
a properly antisymmetric | wave function so that one should be suspicious 
of particle-core coupling models for the | state which include only the lowest 
O"*" and 2"*" core states. In fact, the shell- model parentage is considerably larger 
to the second 2+ state than to the first, as one might guess from the Nilsson 
model picture of the states involved; the spectroscopic factors for the Oj^, 2]*", 
and 2t states of the core are 0.75, 0.45 and 1.00 (0.63, 0.65, and 0.93 for the 
(8-16)2BME interaction of Cohen and Kurath [38]). In contrast the parentage 
for the 1^ ground state of ^^Be is mainly to the first 0+ and 2+ states with 

i^L = 0. 

That many of the known states of ^^Be are mainly of {sdY character was 
suggested long ago and is rather obvious from an analysis of ^''Be(t,p) data 
and the resultant comparison with known (sd)^ states in ^^C and ^^C [66]. For 
^^Be and ^''Be, neutrons have to occupy sd orbits. 



4 Non-normal pEirity states in the p-shell 

States with a nominal Ifko excitation energy are formed by Os — > Op and Op — > 
IsOd excitations with both configurations usually being present to prevent 
unphysical centre-of-mass excitations but with the latter dominating at low 
excitation energies. One interaction that has been widely used one due to 
Millener-Kurath (MK) developed for a study of the j3 decay of ""^^B [67]. The 
MK interaction was "hand crafted" to give a good account of the relative 
separation of the lsi/2 and 0^5/2 centroids as a function of mass number, 
through ^^C and down to ^^Be. More attention was given to ^4 = 11 by Teeters 
and Kurath [68] and by Millcncr and collaborators in studies of the the (3 decay 
of ^^Be [69] and the fast El transition in -"^^Bc [70]. Warburton and Brown [24] 
have subsequently shown that it is possible to obtain a successful fit with an 
rms deviation of ~ 330 keV to a large number of cross-shell energies from 
^4 = 10 — 22, including some for 2hL0 configurations. 

The lowering of the lsi/2 orbit with respect to the 0^5/2 orbit is expected for 
a simple potential well. In shell-model calculations it arises because in the 
potential energy contributions to the single-particle energies at ^^O, the 0^5/2 
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nucleoli interacts more strongly with the Op nucleons while the lsi/2 miclcon 
interacts more strongly with the Os nucleons. The sum of the interactions with 
both shells of the core is similar at ^ "^O but the Od orbit loses attraction relative 
to the Is orbit as p-shell nucleons are removed. In fact, half the difference 
between the interactions with the full p shell roughly accounts for the 2.65 
McV shift in the relative I"*", |^ separation between ^^O and ^^Be. The shift 
between ^''O and ^^C is related simply to the properties of the T=l 
and interactions which may be read directly from the spectrum of ^^N 

[10]. These four particle-hole matrix elements control the basic energetics of 
all the heavy carbon and nitrogen isotopes (see Section 5). 

Table 2 shows the systematics of energies with respect to neutron thresholds, 
"Coulomb energies" , decay widths, and dominant weak-coupling components 
of the proton-rich members of isospin multiplets for Ihw states from A ~ 
11 — 17 (one exception is that for A — IA the widths are given for ^^N because 
of the paucity of information, apart from the 1~ state, for ^^O). All the states 
have isospin T = Tc + |; states with lower isospin gain considerable binding 
energy because a higher spatial symmetry is possible. Because the sd-shell 
nucleon is not restricted by the Pauli Principle, the neutron separation energies 
show a rather smooth variation with mass number (third column of Table 
Table 2), in analogy with the classic example of A separation energies [71]. The 
separation energy eventually has to cross zero as the mass number decreases. 
At the same time, the separation energy for a p-shell neutron is going towards 
zero as neutrons are added so that there is inevitably competition between (at 
least) Qhu, Ihio, and 2fvjj configurations. The energies of the lowest [A — 2)® 
i/(sd)^(0+) states relative to the two-neutron separation threshold also show 
a very smooth behaviour, with approximately constant values along isotopic 
chains; e.g., (6.53, 6.69, 5.47) MeV for ^^-i^C, (4.54, 4.57, 3.74) MeV for 
and (2.30, 3.91, 3.67) MeV for lo-^^gg 

The Coulomb energies given in Table 2 are not true Coulomb energy differences 

except for T = ^. Rather they reflect the differences between neutron and pro- 
ton separation energies at cither end of the isospin multiplet. The calculated 
AEq for a pure weak-coupling configuration involves choosing the well depth 
for a Woods-Saxon well of standard geometry to fit the neutron separation 
energy, then turning on the Coulomb potential of a uniformly charged sphere 
and computing the proton separation for N and Z interchanged. For mixed 
configurations, the calculated AEq- are weighted by the shell-model parent- 
ages, enabling the binding energy of the proton-rich member to be estimated. 

To take the example of the ^^Be, ^^N pair, AE^ is 1.58 MeV for a pure lsi/2 
configuration and 2.92 MeV for a 2f (E) configuration (15% of the shell- 
model wave function). Ai?^ is slightly higher still for the remaining 4% of the 
shell-model wave function. Performing the average puts the ^^N ground state 
1.34 MeV above the ^°C -\- p threshold. This can be compared with a recent 
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Table 2 

Observed and calculated Coulomb energy shifts and nucleon decay widths for IHu 
states which have dominantly lsi/2 or 0(1^/2 character. The trend with changes in 



binding 


enerj 


5y is made clearer by normalizinj 


g all A£ 


'C" to a Z> 


= 6 core. 


Nucleus 




En 






Yth 




Comments 


iiBe 


1 + 

2 


-0.503 


1.845 


1.773 


1467 


1440(200) 


81% si/2{gs) 


12b 


2- 


-1.696 


2.309 


2.286 


87 


118(14) 


70% si/2(5s) 




1- 


-0.749 


1.942 


1.948 


894 


750(250) 


76% si/2{gs) 




1- 


-1.194 


2.046 


~ 2.09 


311 


260(30) 


77% 5i/2(r) 


13c 


1 + 

2 


-1.857 


2.238 


2.278 


32 


31.7(8) 


89% si/2(<7s) 




3 + 
2 


-1.699 


2.189 


2.205 


84 


115(5) 


87% Si/2(2+) 




5 + 
22 


-2.521 


2.437 


2.502 


7 


1 1 


(4/0 Sl/2l,^ } 


14c 


1" 


-2.083 


2.280 


2.252 


28.8 


30(1) 


76% si/2(5s) 


15c 


1 + 
2 


—1.218 


1.857 


2.016 


934 


1000(200) 


98% si/2(5s) 


16n 


0- 


-2.371 


2.168 


2.180 


21 


40(20) 


100% si/2(5s) 


16n 


1- 


—2.094 


2.113 


2.117 


78 


< 40 


97% si/2(5s) 


1^0 


1+ 

2 


-3.273 


2.301 


2.376 






100% si/2(<7s) 


iiBe 


5 + 
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1.275 
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2.475 
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67% d^/2{gs) 
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70% 4/2(1") 


13c 


9 + 
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91% 4/2(2+) 




7 + 
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-1.892 


2.731 


2.666 


3 


9.0(5) 


92% 4/2(2+) 




5 + 
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-1.093 


2.675 


2.695 


45 


47(7) 


80% 4/2 (5s) 
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3- 
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2.636 


2.651 
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3.5 
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< 15 


96% 4/2 (5s) 
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100% 4/2(55) 
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result of 1.27l|j:J^ MeV from resonance elastic scattering of protons [72] (see 
also a result from transfer reactions [73]). Although the complex resonance 
energy is close to becoming unstable, the computed width in Table 2 is close 
to the value of 1.44(20) MeV from the same experiment. All in all, Table 2 
shows remarkable agreement with experiment for both energies and widths 
even in the case of broad states. Encouraged by the similarity in structure 
and separation energy for -^^N(lj") and -'^^N(|'''), one could also extrapolate 
the experimental enegies in Table 2 for low lsi/2 neutron binding energy to 
make a reliable estimate for the ^^N ground state energy. 

The differences in Coulomb energies for the Op, Is, and Od orbits, which are 
considerable, and their behaviour as a function of binding provide a very sen- 
sitive tests of the structure of nuclear states. The well-known Nolan-Schiffer 
anomaly for Coulomb energies can be bypassed by taking only direct Coulomb 
energies [74] with an overall energy scale set by the radius parameter of the 
single-particle well, as has been done in Table 2. The agreement with experi- 
ment suggests that the shell-model parentages are realistic, including the case 
of ^^Be over which there has been much debate. A recent measurement of 
II — — 1.6816(8)/iAr for the magnetic moment of ^-"^Be [75], which is sensitive 
to the 2+ (g) o?5/2 admixture [76], tends to confirm the MK parentage which 
gives fi = — 1.71/XAr when bare-nucleon g-factors are used. Note that Suzuki 
et al. [76] used a a quenched spin g-factor based on the tabulated value of 
H = — 1.315(70)/iAr for ^^C. However, this value appears only in a conference 
proceedings and I strongly suspect that it is not correct; based on the parent- 
age for ^^C ground states in Table 2, //(^^C) should be closer to the Schmidt 
value than /x(^^Be). 

Recent experiments come to the same conclusion about the parentage of the 
^^Be ground state. Analysis of the one-neutron knockout reaction with coinci- 
dent detection of 7 rays from the ^°Be core [77] determines that the ground- 
state spectroscopic factor of 0.74 from the WBP interaction [24] is consistent 
with the data; the slightly higher value from the MK interaction would give 
even better agreement. Likewise, analysis of the ^^Be{p,dY^Be reaction in 
inverse kinematics leads to a 2+ (g) ^5/2 admixture of ~ 16% [78]. 

In the one-nucleon knockout reaction (^^Be,^^ Be + 7) reaction on a '^Be target 
[23], comparable spectroscopic factors are extracted for the bound | and 

I states of ^^Be explicitly demonstrating the breakdown of the = 8 shell 
closure. The shell-model calculations using the WBP interaction [24] reproduce 
the data quite well. They also predict substantial parentage to the unbound 
1^ state of ^^Be [23]. The latter parentage is a warning of the inadequacy, over 
and above center of mass problems, of calculations which limit the number of 
sd orbits to the lsi/2 orbit. Even the ^3^2 configuration is important because 
of a large matrix element with 0^5/2 • Proton removal from ^^B would provide 
another way to study the mixing of Ohu; and 2fiu; 0+ and 2+ states in ^^Be, 
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the predicted C'^S values for pickup to the lowest p-shell 0+ and 2"'" being 
0.52 and 2.12. respectively. Also in -"^^Be, a 1^ state has been seen at 2.68(3) 
MeV by observing 7 rays following the inelastic scattering of ^^Be on a lead 
target [79]. 

The parentage of ^^Be and ^^Be to ^^Be is clearly complicated by the fact that 
{sd)"^ configurations form a substantial component of low-lying states in ^^Be. 
A state in ^^Be with probable 0^5/2 parentage to ^^Be has been located 2.0 
MeV above the neutron threshold [80]; hoth. p^{sd) and [sdf' configurations 
probably need to be considered. A number of other states have also been 
seen [80,81] including a possible ground state 0.78 MeV above threshold with 
^^Be(|^ ) ® {sdy as the most likely candidate. 



5 The heavy boron, carbon, and nitrogen isotopes 

From the discussion in the previous section, it is clear that the "stretched 
isospin" wave functions for a p-shell core plus neutrons in the sd shell can 
usually be described in terms of a few dominant weak-coupling components 
involving low-lying states of the core. The same is true for cases involving 
p-shell proton holes and sd-shell neutrons because then only the weak and 
repulsive T — 1 particle-hole interaction acts. For the nitrogen and carbon 
isotopes one can assume to a first approximation that the lowest states are 
formed by coupling the ^'^N or ^"^C ground state to states of the appropriate 
oxygen isotope. For p~'^{sdY configurations, an estimate of the mass excess 
can be otained from 

ME{A) = ME{16 - m) + ME{16 + n) - ME{^^0) + {V^h) , (4) 

where the ME denote experimental mass excesses and (Vph) is a weighted 
average of mn particle-hole interactions, dominated by the Pi/^2'^5/2 and Pi^2-^i/2 
interactions, which can be taken from ^^N. For a simple monopole interaction 
of the form a + btp.th, {Vph) = mn{a + 6/4) and the fact that an essentially 
constant value for a + b/A can be extracted from any of the known heavy C 
or N isotopes demonstrates the basic consistency of the approximation. Of 
course there is mixing in a shell-model calculation, a good example being the 
importance of small p^^^ admixtures for /x(^^N(l~)), for which two discordant 
results have been reported [82] . 

Two recent experimental studies involving one-neutron removal reactions have 
used wave functions based on the WBP interaction [24] in their analysis. The 
first is a systematic study on 23 neutron-rich psd-shell nuclei with Z = 5 — 9 
and A — 12 — 25 [83] and the second, with coincident 7-ray detection, is a 
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study of i6,i7,i9q |-g^j latter work, the ground-state spins of ^^C and 

^^C are assigned as |^ and respectively. This means that the particle-hole 
interaction inverts the order of the ^^O ground-state doublet in ^''C but does 
not lower the configuration involving the 1.472-MeV state of ^^O enough 
to become the ground state of ^^C. However, the converse is true in ^^C for 
the configuration involving the 1.33-MeV state of ^^O. 

In the case of the boron isotopes, particle-hole interactions involving a P3/2 
hole necessarily play an important role. The tabulations [13] show no definite 
spin-parity assignments for excited states of the core nucleus ^^B. Recently, 

Q _|_ r -j- 

spins of I and | have been suggested for the 3.48 and 3.68 MeV states on the 
basis of a one-neutron knockout experiment [85]. In addition, a wealth of new 
information has been obtained on states of ^^-isg ^sing multinucleon transfer 
reactions [86]. For example, the strongest state in one- proton removal, two- 
neutron addition reactions on ^^C leading to ^^B is a state at 6.42 MeV. This 
state almost certainly corresponds to the I"*" T=| state at 21.47 MeV in -'■^C 
which is strongly excited by an M4 transition in inelastic electron scattering. 
The predicted width for the ^^B(2"'") ® 1^0^5/2 analog in ^^B is 33 keV, which is 
in excellent agreement with the tabulated width of 36(5) keV. Measurements 
have also been made for the ground-state moments of a number of the boron 
isotopes and successfully interpreted in terms of shell- model calculations [87] . 



6 Discussion 

The structure of neutron-rich light nuclei has been discussed in terms of shell- 
model calculations that include all configurations at a given nfuj level of exci- 
tation. Other models which can better take into account the radial degrees of 
freedom necessary to describe loose binding and clustering effects tend to be 
applied to specific cases rather than to describe the broad sytematics of nuclear 
spectra. In any case, the results of these calculations are usually interpreted 
in terms of the essential shell-model configurations involved. 

The major question that has to be addressed in any calculation is how to treat 
the mixing of configurations which differ in energy by 2hw or more. Certainly, 
if only the lowest few states of each diagonalized nfuv space are kept and mixed 
in what might be called an extended coexistence model, good spectra and a 
good description of the relationships (parentage) between low-lying states in 
neighbouring nuclei can be achieved. This has been the approach in a number 
of the cluster-model calculations referenced earlier. 

The shell-model runs into consistency problems if diagonalizations are per- 
formed in complete {0 + 2)}tio spaces rather than using the "coexistence" ap- 
proach [88,89,24]. The crux of the matter is that the dominant Ahw — 2 
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interaction transforms as a (2 0) SU(3) tensor (just as most of H in Eq. (1) 
transforms as (0 0) for Ahuj = 0) leading to large matrix elements between 
quite widely separated configurations and slow convergence as a function of 
hjjj. Powerful mathematical techniques exist for the symplectic shell model, a 
natural extension of Elliott's SU(3) model, and the cluster model to handle 
calculations in "vertically" truncated model spaces (see articles in Ref. [40]). 
The correlations induced in low-energy wave fucntions can be viewed both 
as improvements to the radial wave functions between clusters and as an in- 
troduction of the RPA-type correlations necessary to satisfy energy-weighted 
sum rules. 
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